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The  entropy  changes  (AS)  in  various  cathode  and  anode  materials,  as  well  as  complete  Li-ion  batter¬ 
ies,  were  measured  using  an  electrochemical  thermodynamic  measurement  system  (ETMS).  A  thermal 
model  based  on  the  fundamental  properties  of  individual  electrodes  was  used  to  obtain  transient  and 
equilibrium  temperature  distributions  of  Li-ion  batteries.  The  results  from  theoretical  simulations  were 
compared  with  results  obtained  in  experimental  measurements.  We  found  that  the  detailed  shape  of  the 
entropy  curves  strongly  depends  on  the  manufacturer  of  the  materials  even  for  the  same  nominal  com¬ 
positions.  LiCo02  has  a  much  larger  entropy  change  than  LiNixCoyMnz02.  This  means  that  LiNixCoyMnz02 
is  much  more  thermodynamically  stable  than  LiCo02).  The  temperatures  around  the  positive  terminal  of 
a  prismatic  battery  are  consistently  higher  than  those  at  the  negative  terminal,  due  to  differences  in  the 
thermal  conductivities  of  the  different  terminal  connectors.  When  all  other  simulation  parameters  are  the 
same,  simulations  that  use  a  battery-averaged  entropy  tend  to  overestimate  the  predicted  temperatures 
when  compared  with  simulations  that  use  individual  entropies  for  the  anode  and  the  cathode,  due  to 
computational  averaging. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High-power  secondary  batteries  for  plug-in  hybrid  electric  vehi¬ 
cle  applications  generate  much  more  heat  during  rapid  charge  and 
discharge  cycles  at  high  current  levels,  such  as  during  quick  accel¬ 
eration  when  the  battery  consumes  power  or  during  deceleration 
or  quick  stops  when  the  battery  is  charged  using  recovered  inertial 
energy  [1  ].  However,  most  previously  used  lithium  batteries  oper¬ 
ated  at  low  discharge  rates.  The  use  of  these  batteries  at  high  power 
rates  will  pose  serious  thermal-stability  problems,  especially  when 
the  battery  size  increases  and  many  batteries  are  interconnected  to 
form  large  packages.  Toyota  recently  announced  that  because  of 
these  safety  concerns,  it  will  delay  by  1-2  years  the  launch  of  new 
high-mileage  hybrid  vehicles  equipped  with  Li-ion  battery  tech¬ 
nology  [2].  As  battery  size  increases,  the  ratio  of  cooling  area  to 
heat  generating  volume  decreases,  and  as  charge/discharge  current 
increases,  more  heat  is  generated.  The  temperature  of  the  battery 
thus  rises  dramatically,  leading  to  the  possibility  that  internal  tem¬ 
peratures  will  exceed  permissible  levels.  The  transient  variation 
of  the  temperature  distribution  in  the  battery  during  operation 
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should  therefore  be  considered  during  development  of  large  Li-ion 
secondary  batteries. 

In  a  model  established  two  decades  ago,  Bernardi  et  al.  [3] 
described  a  general  energy  balance  equation  for  battery  systems. 
This  rather  complete  system  of  equations  laid  the  foundations  for 
the  prediction  of  cell  temperatures  and  heat  generation  result¬ 
ing  from  electrochemical  reactions,  phase  changes,  mixing  effects, 
and  Joule  heating.  In  1993,  Doyle  et  al.  [4]  developed  a  micro¬ 
scale  continuum  model  of  the  galvanostatic  charge  and  discharge 
of  Li/polymer/porous  cathode  cells  using  concentrated  solution 
theory.  In  that  model,  they  simulated  insertion  of  Li  into  the  cath¬ 
ode  using  superposition  of  concentration  fields  to  simplify  the 
calculations.  Concentration  and  charge  distributions  across  the 
anode/electrolyte/cathode  assembly  were  predicted,  but  thermal 
performance  was  not.  Their  work  was  followed  by  a  series  of  papers 
focused  on  thermal  performance.  Assuming  uniform  heat  gener¬ 
ation  throughout  the  battery,  superposition  methods  were  again 
employed  by  Newman  and  Tiedemann  [5]  to  compute  the  tem¬ 
perature  rise  for  a  range  of  aspect  ratios  in  a  three-dimensional 
battery.  Methods  for  treating  time-dependent  heat  generation  rates 
were  also  discussed,  but  material  temperature  dependencies  were 
not  treated.  In  a  later  report,  Pals  and  Newman  [6]  established 
a  one-dimensional  model  for  a  single  cell  and  for  stack  thermal 
performance.  Material  properties  were  temperature  dependent. 
The  entropic  term  in  the  heat  generation  rate  was  recognized,  but 
deleted  from  the  analysis  because  insufficient  data  were  available 
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for  the  computations.  Consequently,  the  heat  generation  rate  was 
essentially  constant  for  each  simulation.  Spatial  distributions  of 
temperature  and  ion  concentrations  were  produced  for  the  stack 
model.  In  1996,  Doyle  et  al.  [7]  returned  to  more  detailed  simula¬ 
tions  of  local  electrochemistry  and  ion  concentration  distributions. 
Others  addressed  the  thermal  performance  of  lithium  ion  batter¬ 
ies.  Funahashi  et  al.  [  8  ]  solved  for  transient  temperatures  using  bulk 
properties,  and  uniform  heat  generation  over  the  battery,  but  only 
Cp  was  temperature  dependent.  Thomas  and  Newman  [9]  found 
the  effects  of  the  heat  of  mixing  to  be  negligible  in  a  well-designed 
battery,  using  a  very  complete,  one-dimensional  analysis.  Araki  et 
al.  [10]  used  a  one-dimensional  model  with  lumped  parameters  for 
NiMH  battery  thermal  predictions.  Onda  et  al.  [  11  ]  modelled  a  one¬ 
dimensional  (radial)  battery  transient  with  lumped  parameters. 
Inui  et  al.  [12]  developed  transient  two-dimensional  and  three- 
dimensional  models  for  temperature,  charge  state,  and  current  as  a 
function  of  position.  Although  their  method  appears  to  have  been  a 
numerical  solution  of  simultaneous  equations,  they  do  not  describe 
the  mathematical  technique  they  used. 

Theoretical  predictions  of  important  physical  properties  have 
also  recently  become  available.  Kim  and  Pyun  [13]  used  a  Monte 
Carlo  approach  to  determine  basic  thermodynamic  quantities 
(internal  energy,  entropy)  for  Li  intercalation  in  Li!_xMn204  elec¬ 
trodes,  obtaining  good  agreement  with  experimental  results. 
Thermodynamic  properties  also  are  available  for  LixCo02  [14] 
and  for  Li  intercalation  into  carbons  [15].  Extensions  to  three- 
dimensional  thermal  performance  of  batteries  have  been  exten¬ 
sively  addressed  by  Pesaran  et  al.  [16,17]  at  the  National  Renewal 
Energy  Laboratory.  They  successfully  simulated  thermal  behav¬ 
ior  of  a  cylindrical  Li-ion  battery  assuming  that  the  core  material 
(cathode/separator/anode)  consisted  of  a  homogenous  material 
with  average  properties  for  resistivity,  thermal  conductivity,  and 
entropy.  They  also  investigated  various  ventilation  methods  for 
battery  cooling.  They  used  the  finite  element  method  to  model 
the  complete  battery  system,  including  external  heat  transfer.  In 
these  studies,  entropy  changes  in  battery  electrodes  were  rec¬ 
ognized  as  an  important  factor  that  may  significantly  affect  the 
thermal  behavior  of  Li-ion  batteries.  For  example,  Hong  et  al. 

[18]  indicated  that  the  entropy  change  may  contribute  more  than 
50%  of  the  total  heat  at  the  C/1  discharge  rate.  However,  the 
entropy  changes  in  previous  studies  have  been  assumed  to  be  uni¬ 
form  parameters  for  complete  batteries,  partly  because  entropy 
data  for  individual  electrodes  were  not  available.  Recently,  Yazami 

[19]  developed  an  approach  for  the  measurement  of  entropy  and 
enthalpy  changes  in  electrode  materials  (such  as  LiCo02,  LiFeP04, 
and  graphite). 

To  predict  more  detailed  thermal  behavior  for  a  battery  with 
different  cathode  and  anode  materials,  it  is  necessary  to  develop 
a  versatile  model  based  on  the  individual  properties  (including 
entropy  changes)  of  the  cathode  and  anode  separately.  This  model 
can  be  used  not  only  to  simulate  the  thermal  properties  of  existing 
batteries,  but  also  to  predict  the  internal  temperature  of  the  bat¬ 
tery  as  a  function  of  battery  chemistry,  state  of  charge  (SOC),  power, 
and  battery  geometries.  Such  an  improved  model  will  help  devel¬ 
opers  in  industry  design  a  new-generation  battery  that  minimizes 
heat  generation  and  thermal  fluctuations,  therefore  significantly 
reducing  heat-induced  failures. 

2.  Simulation  of  the  thermodynamic  behavior  of  lithium 
ion  batteries 

2.1.  Model  description 

The  thermodynamic  properties  (entropy  change,  AS)  of  cath¬ 
odes  and  anodes  materials  published  recently  by  Yazami  [19]  were 


used  in  an  initial  effort  to  simulate  the  thermal  behavior  of  bat¬ 
teries.  The  transient  temperature  distribution  of  the  battery  can  be 
solved  from  the  following  differential  equation: 

p(x)Cp(x,  T(x,  t =  — V[— k(x,  T(x,  t))VT(x,  t)]+Q(x,  T(x,  t)) 

(1) 


where 

Q(x,  T(x,  t))  =  R(x,  T(x,  t))/(x,  t)2  -  T(x,  t)AS(x,  T(x,  t))^P  (2) 

nF 

is  the  space-  and  temperature-dependent  heat  generation  rate.  The 
model  variables  are  time  dependent  through  their  dependence  on 
the  temperature  and  current,  which  are  time  dependent.  It  should 
be  noted  that  the  effects  of  the  heat  of  mixing,  heat  effects  due  to 
possible  chemical  reactions,  and  convective  effects  are  neglected 
in  the  model.  The  first  term  in  Eq.  (2)  is  the  rate  of  heat  generated 
from  ohmic  heating.  The  second  term  is  the  rate  of  heat  generated 
from  micro-structural  (entropic)  changes.  The  rate  of  heat  gener¬ 
ated  in  the  electrode  or  battery  can  be  obtained  by  integrating  Eq. 
(2)  over  the  electrode  or  battery  volume.  To  simplify  the  compari¬ 
son  with  experimental  results  (which  normally  plots  temperature 
of  the  battery  as  a  function  of  the  SOC),  in  the  simulations  we  used 
depth  of  discharge  (DOD)  (DOD  =  1  -  SOC)  as  the  independent  vari¬ 
able  to  replace  time  t  in  Eq.  ( 1 ).  We  took  this  approach  because  DOD 
varies  linearly  with  time  (t)  under  constant  the  current  conditions 
used  in  our  experiment  (DOD  =  — Jt/C),  where  C  is  the  battery  capac¬ 
ity  (As)  and  J  is  the  constant  current  (A).  The  differential  equation 
employed  in  the  simulations  was 


J  8T 
CdDOD 


J_d_ 
rP  dr 


dT 

dr 


dT 
ai  dz 


+  Q(T,  DOD) 


ki 


0) 


where  i  denotes  the  material  (package,  anode,  cathode,  etc.); 
a  =  k/pCp  is  the  thermal  diffusivity  for  each  material;  =  0  or  1  for 
rectangular  or  cylindrical  coordinates,  respectively;  z  is  the  axial 
coordinate;  r  is  the  thickness  coordinate;  T  is  temperature;  and  Q 
is  the  heat  generation  rate  for  each  material: 

Qi  =  I2Ri-TAsT  (4) 

where  I  is  the  current  density  (A cm-2),  Rz  is  the  material  resistivity 
(Q-cm),  T is  temperature  (K),  A5  is  entropy  change  (J  mole-1  K-1 ), 
n  equals  one  electron  per  reaction,  and  F  is  the  Faraday  constant 
(96,485  C  mole-1 ).  Finite  difference  methods  for  solving  Eq.  (3)  have 
been  investigated  by  Melgaard  and  Sincovec  [20],  including  the 
body,  edges,  and  corners  of  the  battery.  It  is  perhaps  noteworthy 
that  the  heat-transfer  equation  they  treated  accounts  for  spatial 
variations  in  the  thermal  diffusivities,  a.  The  above  equations  can 
be  applied  to  the  cathode,  separator,  and  anode,  respectively,  and 
matched  by  boundary  conditions.  In  a  typical  LiCo02/C  battery,  the 
overall  electrochemical  reaction  can  be  expressed  as 


2Lio.5Co02  +  LiCg  2Lio.5+xCo02  +  Lii_2xCg  (5) 

The  corresponding  entropy  change  in  this  battery  can  be 
expressed  as 

ASceii  =  ASLio5+xco02  +  ASLi,_2xC6  (6) 

The  total  heat  generated  from  the  LiCo02  cathode  and  the  car¬ 
bon  anode  from  their  micro-structural  changes  can  be  calculated 
from  Eqs.  (4)  and  (6).  Fig.  1  shows  the  entropy  changes  in  a  typical 
LiCo02/C  battery  in  various  charge  states  (entropy  data  are  taken 
from  Refs.  [14]  and  [19],  respectively).  The  term  x  on  the  x-axis 
of  Fig.  1  represents  the  electrode-averaged  ion  concentration  in 
LiCo02  and  carbon  electrodes.  It  is  important  to  notice  that  the 
spikes  in  the  entropy  changes,  which  correspond  to  the  spikes  in 
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Fig.  1.  Entropy  changes  in  a  typical  UC0O2/C  battery  in  various  charge  states. 
Entropy  data  for  (a)  LiCoCh  and  (b)  graphite  are  taken  from  Refs.  [13]  and  [18], 
respectively. 


heat  generation,  in  the  cathode  and  anode  strongly  depend  on  the 
structural  properties  of  the  materials  [9].  Apparently,  if  both  ther¬ 
mal  spikes  in  the  cathode  and  anode  in  a  battery  appear  at  the 
same  time,  the  battery  will  experience  large  thermal  stresses  that 
can  lead  to  early  failure  or  explosions.  However,  appropriate  selec¬ 
tion  of  anode/cathode  pairs  and  their  molar  ratios  can  effectively 
avoid  the  simultaneous  appearance  of  anode  and  cathode  spikes  in 
a  battery,  therefore  eliminating  an  important  source  of  overheat¬ 
ing  during  battery  operation.  The  objective  of  this  paper  is  to  enable 
the  ‘virtual  simulation’  of  such  selections  by  computer  simulations, 
thus  avoiding  the  expense  of  actually  building  many  new  battery 
prototypes  during  material  development  efforts. 

2.2.  Boundary  conditions 

In  Cartesian  coordinates  (/3  =  0),  the  boundary  conditions  can  be 
expressed  generally  as 

„r  +  /,f.c  (7) 

where  x  equals  r  or  z  in  the  thickness  direction,  and  a,  b ,  and  c  are 
at  least  piecewise  continuous  functions  of  the  spatial  dimensions. 
Fig.  2  shows  schematic  of  boundary  conditions  for  an  x-z  section 


Fig.  2.  Schematic  of  boundary  conditions  for  an  x-z  section  through  a  prismatic 
battery.  The  large  flat  face  of  the  battery  is  normal  to  the  plane  of  the  page. 


through  a  prismatic  battery  (not  to  scale).  The  large  flat  face  of  the 
battery  is  normal  to  the  plane  of  the  page.  Eq.  (7)  is  a  general  form 
that  encompasses  all  three  types  of  boundary  conditions.  Exam¬ 
ples  of  straightforward  applications  are:  when  a  =  0  the  boundary 
conditions  are  von  Neumann  type;  when  b  =  0,  the  boundary  condi¬ 
tions  are  Dirichlet  type;  and  when  a  and  b  are  both  nonzero,  Cauchy 
conditions  prevail.  For  a  full  two-dimensional  model  of  the  axial 
and  thickness  directions  for  a  prismatic  battery,  including  the  elec¬ 
trical  leads,  the  following  boundary  conditions  are  used.  Natural 
convective  cooling  by  air  is  assumed  to  dominate  cooling  on  the 
exterior  surface  of  the  battery.  For  a  convective  heat-transfer  coef¬ 
ficient,  hcool,  and  an  ambient  air  temperature,  Tair,  the  boundary 
parameters  are  a  =  -hcool,  b  =  -kri ,  c=-hC00lTair  at  the  right  hand 
battery  face,  where  kri  is  the  thermal  conductivity  of  the  package 
in  the  radial  (thickness)  direction.  At  the  left  face,  the  sign  of  b  is 
reversed  for  mathematical  self-consistency  of  the  temperature  gra¬ 
dients.  The  convective  heat-transfer  coefficient  for  laminar  flow  of 
air  along  a  vertical  surface  is  [21  ]: 

hcool  =  c(~r)  (8) 

where  C  and  n  are  constants,  L  is  the  height  of  the  battery,  and  the 
temperature  difference  is  again  from  the  surface  to  the  ambient  air. 
The  effects  of  the  battery  electrical  leads  are  also  included  on  the 
upper  edge  of  the  battery,  by  conductive  heat-transfer  pathways 
located  just  inside  the  battery  package  layers.  The  pathway  for  the 
anode  is  copper,  while  the  pathway  for  the  cathode  is  aluminum. 
In  this  case,  a  =  -/<M/Ax,  b  =  -/<ze,  and  c=-kM  TWire/Ax ,  where  kM  is 
the  axial  conductivity  of  the  appropriate  lead  wire,  and  Twire  is  the 
measured  temperature  of  the  lead  wire  at  a  distance  of  Ax  from  the 
battery  edge. 


2.3.  Material  properties 


The  material  properties  (e.g.,  thermal  conductivities,  electri¬ 
cal  resistivities,  etc.)  differ  greatly  between  the  various  materials, 
which  can  lead  to  numerically  ‘stiff  problems  and  slow  com¬ 
putations.  As  in  most  numerical  simulations  of  multi-layered 
lithium  batteries,  the  properties  are  ‘lumped’  into  electrode 
units  with  effective  properties  that  exhibit  less  variation.  In  the 
present  case,  the  repeating  electrode  units  consisted  of  separa¬ 
tor/electrode/metal/electrode  assemblies.  The  effective  properties, 
such  as  thermal  diffusivities,  were  calculated  by  volume  averages 
for  both  the  thickness  (i.e.,  the  direction  of  the  current  flow)  and  the 
axial  directions.  Additional  details  are  given  below.  Effective  electri¬ 
cal  resistivities  in  the  radial  (conduction)  direction  were  calculated 
for  electrode  assemblies  (i.e.,  separator/electrode/metal/electrode), 
assuming  that  the  porous  separator  was  dominated  by  the  elec¬ 
trolyte  resistivity,  using  the  following  equation: 


^radial 
'  assembly 


Er<f< 


(9) 


where  q  is  the  thickness  of  the  layer  in  the  assembly,  rz  is  the  elec¬ 
trical  resistivity  of  the  solid  material  for  the  electrode  and  metal, 
and 


r  _  ^electrolyte  rsolid 

1 ~  p  + 1 -p  1 

for  the  separator,  where  p  is  the  porosity.  Electrical  resistivity  in  the 
axial  (z)  direction  is  dominated  by  the  metal  and  is  very  low  com¬ 
pared  to  the  other  materials.  An  analogous  ‘conduction  in  series’ 
equation  is  used  for  the  electrode  assembly  thermal  conductivity 
in  the  radial  direction,  while  the  thermal  conductivity  (/<)  in  the 
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Table  1 

The  effective  material  properties  in  a  lithium  ion  battery. 


axial  direction  (conduction  in  parallel)  is 


/.axial 

'^assembly 


Eli 

E  f.A 


(u) 


Assembly  densities  (p)  are  calculated  as  for  the  radial  assembly 
resistivities,  while  the  specific  heats  (Cp)  are  calculated  by  mass 
averages: 


_  SfiP/cp,i 

Cp, assembly-  ^  f.  p. 


(12) 


The  effective  properties  for  prismatic  lithium  ion  batteries  pro¬ 
vided  by  Manufacturer  A  are  listed  in  Table  1.  “Wet”  properties 
represent  the  material  properties  after  they  are  wetted  with  elec¬ 
trolyte. 


3.  Experiment 

3  A.  Entropy  measurement 

Entropy  data  of  cathodes,  anodes,  and  complete  cells  were 
obtained  by  monitoring  the  change  in  voltage  as  a  function  of  tem¬ 
perature  using  the  ETMS  (Viaspace,  ETMS-1000).  Several  cathode 
and  anode  powders  were  obtained  from  LICO  Technology  Corp. 
Cathodes  and  anodes  were  mixed  with  85%  active  material,  10% 
PVDF,  and  5%  carbon  black,  and  then  cast  on  Al  and  Cu  foil,  respec¬ 
tively.  Customer-coated  cathode  and  anode  films  obtained  from 
Manufacturer  A  were  used  as  received  for  entropy  measurement.  To 
measure  entropy  changes  in  cathodes  or  anodes,  they  were  assem¬ 
bled  in  a  Type-2325  coin  cell  with  a  0.75-mm-thick,  Li-foil  anode 
(Aldrich),  a  Celgard  2501  polypropylene  separator  and  electrolyte 
(1  M  LiPFg  in  EC/DMC  [Merck]),  using  a  coin  cell  crimper  (National 
Research  Council,  Canada).  Complete  cells  with  LiCo02-based  cath¬ 
odes  and  graphite-based  anodes  from  Manufacturer  A  also  were 
assembled  using  a  similar  procedure  to  that  described  above  by 
replacing  Li  metal  with  the  graphite-based  anode. 

The  cathode  cells  were  cycled  at  C/4  rate  between  3  and  4.1  V, 
with  the  cell  potential  held  at  the  end  values  for  2  h  to  complete 
the  charge  or  discharge.  The  voltage  range  for  anode  cells  was 
0.02-0.5  V.  After  three  cycles,  the  cells  were  brought  to  a  fully  dis¬ 
charged  state  prior  to  the  start  of  entropy  determination.  Charge 
stabilization  was  done  at  35  °C  for  2h,  followed  by  changing  the 
cell  temperature  to  the  initial  desired  value.  Typically,  this  initial 
temperature  is  set  to  be  equal  to  charge  stabilization  temperature. 
After  15  min  at  this  initial  temperature,  the  cell  voltage  was  mea¬ 
sured,  and  then  the  cell  temperature  was  changed  to  the  next  value 
(e.g.,  25  °C).  After  temperature  equalization  at  that  temperature, 
the  cell  voltage  was  measured  at  the  end  of  the  period,  and  then 
the  temperature  was  changed  to  the  next  level  (e.g.,  15  °C).  The 
cell  temperature  was  then  brought  to  the  initial  value  to  deter¬ 
mine  deviation  from  the  initial  voltage  reading.  Entropy  changes 
were  calculated  from  the  slope  of  the  voltage-time  curve  using 


the  equation  AS  =  F[8AF/8T),  where  F  is  the  Faradaic  constant.  The 
experiments  were  repeated  for  LiCo02 /graphite  cells,  using  voltage 
limits  of  4.2  and  2.5  V. 

3.2.  Temperature  profile  measurement  during  cycling  of  prismatic 
cells 

Prismatic  lithium  ion  batteries  (using  modified  LiCoO 2/ 
graphite)  with  nominal  capacity  of  4.5  Ah  from  Policell  Technolo¬ 
gies,  Inc.  were  used  for  temperature  measurements  under  different 
conditions.  These  batteries  have  eight  cells  stacked  inside  and 
have  a  dimension  of  10  cm  x  10  cm  x  0.33  cm.  These  batteries  were 
cycled  by  charging  at  C/3  rate  to  4.2  V  and  discharging  at  C/3,  C, 
2C,  and  4C  to  2.8  V  using  a  BTS  500  battery  cycler  (Firing  Circuits). 
The  temperature  on  each  face  of  a  battery  was  monitored  using 
type  K  thermocouples.  A  schematic  location  of  the  thermocouples 
is  shown  in  Fig.  3.  The  battery  was  suspended  in  a  Plexiglass  frame 
to  allow  air  exposure  on  all  sides  for  convective  cooling.  Small  alu¬ 
minum  clamps  were  attached  to  its  electrode  leads.  Thermocouples 
were  attached  to  both  faces  of  the  battery  as  shown  in  Fig.  3. 

4.  Experimental  and  simulation  results 

In  our  initial  simulation,  the  entropies  of  LiCoO 2/C  reported  in 
Refs.  [14,19]  have  been  used  to  compare  the  effect  of  the  entropy 
contribution  to  the  battery  temperature  when  (1)  the  average 
entropy  of  the  complete  battery  was  used  and  (2)  the  individual 
entropies  of  the  anode  and  cathode  were  used  in  the  simulations. 
Fig.  4  shows  the  temperature  profile  simulation  for  an  eight-cell 
LiCo02/C  prismatic  battery.  For  clarity,  only  the  temperatures  that 
correspond  to  point  5  in  Fig.  3  have  been  shown  in  Fig.  4.  Tcenter  and 
^surface  represent  the  temperatures  at  the  center  of  eight-cell  stack 
and  on  the  surface  of  the  package.  The  battery  capacity  was  4.8  Ah, 
the  convective  cooling  coefficient  was  0.001  W cm-2  °C_1,  and  the 
ambient  room  air  and  battery  initial  temperature  was  25  °C.  Fig.  3 
indicates  that  the  effect  of  using  of  the  battery-averaged  entropy 
was  to  overestimate  the  predicted  temperatures.  When  all  other 
simulation  parameters  are  the  same,  the  temperature  profile  of  the 
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Fig.  3.  Thermocouple  locations  in  a  prismatic  lithium  ion  battery. 
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Fig.  4.  Temperature  profiles  in  a  UC0O2/C  battery  with  an  eight  bi-cell  stack  at 
various  discharge  states  when  discharged  at  the  1C  rate. 


Fig.  6.  Entropy  changes  for  modified  LiCo02  (vs.  Li)  prepared  by  Manufacturer  A. 
Two  curves  in  the  figure  represent  the  measured  results  on  different  days  for  the 
same  cell. 


battery  is  about  5  °C  higher  (during  the  later  stages  of  the  discharge) 
when  the  average  battery  entropy  was  used  when  compared  to  the 
profile  obtained  when  separate  entropies  are  used  for  the  anode 
and  the  cathode.  The  reason  for  this  difference  is  because  different 
methods  of  averaging  material  properties  result  in  different  mate¬ 
rial  properties.  For  best  accuracy,  it  is  recommended  that  material 
properties  be  averaged  over  volumes  as  small  as  possible.  It  also  is 
interesting  to  note  that  no  significant  temperature  difference  was 
observed  between  the  center  (along  the  battery  thickness  direc¬ 
tion  and  inside  of  the  battery)  of  the  stack  and  on  the  surface  of  the 
battery  when  this  eight-cell  battery  was  discharged  at  the  1C  rate 
under  natural  ventilation  conditions. 

To  obtain  a  better  comparison  between  the  simulation  and 
experimental  data,  entropy  changes  for  cathodes,  anodes  versus 
lithium  metal  and  for  complete  cathode/anode  combinations  have 
been  measured  as  a  function  of  the  SOC  using  the  ETMS  (Viaspace, 
ETMS-1000).  Fig.  5  shows  A5  versus  SOC  for  LiCo02  obtained  from 
LICO  Technology  Corporation.  Fig.  6  shows  AS  versus  SOC  for  mod¬ 
ified  LiCo02  obtained  from  Manufacturer  A.  The  two  curves  in  the 
figure  represent  the  measured  results  on  different  days  for  the  same 
cell.  In  comparing  the  AS  of  LiCo02  shown  in  Fig.  1(a)  with  those  in 
Figs.  4  and  5  (100%  SOC  at  the  right  side  of  Figs.  4  and  5  corresponds 
to  x  =  0.5  at  the  left  side  of  Fig.  1),  it  is  interesting  to  note  that  all 
of  the  LiCo02-based  cathodes  show  a  valley  in  their  entropy  data 
at  10-20%  SOC  or  x  ~  0.9.  However,  no  peak  at  x  =  0.55  or  SOC  =  90% 
was  observed  in  Figs.  5  and  6.  The  difference  between  the  entropy 


Fig.  5.  Entropy  changes  for  commercially  available  LiCo02  cathode  (LICO  Technology 
Corporation). 


curves  obtained  from  the  LiCo02 -based  cathode  may  be  attributed 
to  the  different  manufacturing  processes  and  impurity/doping  used 
by  different  manufactures. 

Fig.  7  shows  AS  versus  SOC  for  LiNixCoyMnz02  obtained 
from  LICO  Technology  Corporation.  The  agreement  between  two 
runs  (shown  as  two  separate  curves  in  the  figure)  on  differ¬ 
ent  days  for  the  same  cell  is  reasonably  good.  It  is  of  interest 
to  notice  that  the  entropy  changes  of  the  LiNixCoyMnz02  cath¬ 
ode  shown  in  Fig.  7  are  largely  different  as  compared  with  the 
LiCo02  cathode  shown  in  Figs.  5  and  6.  First,  the  overall  entropy 
change  in  LiNixCoyMnz02  is  only  6Jmol-1  K-1.  The  correspond¬ 
ing  enthalpy  change  is  43.2 J  mol-1  K_1.  However,  Figs.  5  and  6 
show  that  the  overall  entropy  change  in  LiCo02  is  more  than  80 
and  70  J  mol-1  K-1 ,  respectively.  The  larger  entropy  change  (mainly 
attributed  to  the  change  in  configurational  entropy  [14])  represents 
a  larger  structure  instability  in  LiCo02  than  in  LiNixCoyMnz02.  We 
also  noted  that  the  largest  entropy  change  in  these  two  materi¬ 
als  appears  at  two  different  charge  states.  For  LiNixCoyMnz02,  the 
largest  entropy  change  appears  when  SOC  =  0  (i.e.,  at  the  fully  dis¬ 
charged  state).  For  LiCo02,  the  largest  entropy  change  appears  at 
SOC  ~  20%. 

Fig.  8  shows  AS  versus  DOD  (=1  -  SOC)  for  a  graphite-based 
anode  obtained  from  Manufacturer  A.  The  entropy  data  for  graphite 
shown  in  Figs.  1(b)  and  8  are  very  similar.  In  this  case,  the  x  axes 
on  the  left  sides  of  both  curves  represent  the  fully  lithiated  state  of 


Fig.  7.  Entropy  changes  for  commercially  available  LiNixCoyMnz02  (LICO  Technology 
Corporation).  The  agreement  between  two  runs  (shown  as  two  separate  curves  in 
the  figure)  on  different  days  for  the  same  cell  is  reasonably  good. 
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Fig.  8.  Entropy  changes  for  a  graphite-based  anode  (vs.  Li)  prepared  by  Manufacturer 
A.  The  two  curves  in  the  figure  represent  the  measured  results  on  different  days  for 
the  same  cell. 


graphite.  The  two  curves  in  Fig.  7  represent  the  measured  results  on 
different  days  for  the  same  cell.  The  entropy  data  for  the  complete 
coin  cell  using  modified  LiCo02  and  graphite-based  anode  prepared 
by  Manufacturer  A  is  shown  in  Fig.  9. 

Temperature  profiles  of  prismatic  lithium  ion  batteries  during 
charge  and  discharge  were  measured  for  two  conditions.  Fig.  10 
shows  temperature  profiles  of  the  samples  at  the  center  of  the  sam¬ 
ple  (point  5  in  Fig.  3)  for  an  un-insulated  0.5  Ah  Li-ion  cell  charged 
at  the  C/3  rate,  and  the  discharged  at  C/3,  C,  and  2C  rates.  A  small 
increase  in  temperature  was  observed  at  the  end  of  charge.  This 
increase  probably  was  the  result  of  the  high  impedance  of  the 
battery  at  extreme  states  of  charge.  The  same  pattern  was  evi¬ 
dent  during  discharge  at  extreme  states  of  charge.  Interestingly, 
immediately  after  discharge  when  charging  was  commenced,  the 
temperature  decreased  because  charging  for  the  Li-ion  battery 
is  endothermic.  This  is  consistent  with  the  thermal  behavior  of 
lithium  ion  batteries  reported  by  Flong  et  al.  [18].  Fig.  11  shows  the 
temperature  profiles  at  the  center  of  the  sample  (point  5  in  Fig.  3) 
for  an  insulated  4.5  Ah  Li-ion  cell  charged  at  the  C/3  rate,  and  the 
discharged  at  C/3,  C,  and  2C  rates.  In  comparisons  of  Figs.  10  and  11, 
we  found  that  the  thermally  insulated  samples  had  a  much  slower 
temperature  response,  but  the  maximum  temperature  (45.5  °C)  of 
the  sample  is  higher  than  those  of  the  air-cooled  samples  (36.4  °C). 


Fig.  9.  Entropy  changes  for  modified  LiCo02  vs.  a  graphite-based  anode  prepared  by 
Manufacturer  A  after  cycling  at  the  C/2  rate  between  2.9  and  4.1  V.  The  two  curves 
in  the  figure  represent  the  measured  results  on  different  days. 


Fig.  10.  Temperature  profile  at  the  center  (point  5  in  Fig.  2)  of  an  un-insulated  4.5  Ah 
Li-ion  cell  charged  at  the  C/3  rate,  and  then  discharged  at  the  C/3,  C,  and  2C  rates. 


Fig.  11.  Temperature  profile  of  an  insulated  4.5  Ah  prismatic  Li-ion  cell  charged  at 
the  C/3  rate,  and  then  discharged  at  the  C/3,  C,  and  2C  rates. 

The  surface  temperature  variation  of  an  un-insulated  sample 
(4.5  Ah  prismatic  battery  obtained  from  Policell  Technologies,  Inc.) 
was  measured  with  multiple  thermocouples  (see  Fig.  3  for  the  loca¬ 
tion  of  the  thermocouples).  Fig.  12  shows  the  temperature  profile 
of  the  sample  when  discharged  at  the  3 C  rate.  We  found  that  the 
temperatures  (at  points  1,  2,  10,  11)  around  the  positive  terminal 
(i.e.,  the  Al  tab)  are  consistently  higher  (~5°C)  than  those  at  the 


Time/s 

Fig.  12.  The  surface  temperature  profile  of  an  un-insulated  prismatic  sample 
(4.5  Ah)  (see  Fig.  2  for  the  location  of  the  thermocouples).  The  sample  was  discharged 
at  the  3C  rate. 
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Fig.  13.  Experimental  and  simulated  temperature  profiles  for  an  eight-cell  prismatic 
battery  discharged  at  the  4 C  discharge  rate. 

negative  terminal  (i.e.,  the  Cu  tab).  This  finding  is  consistent  with 
those  observed  by  Pesaran  et  al.  on  cylindrical  cells  when  the  two 
terminals  were  on  the  same  side  [16].  We  also  found  that  the  sam¬ 
ple  temperature  decreased  from  the  top  (close  to  the  tab  area)  to 
the  bottom  of  the  sample.  Minimal  temperature  variations  were 
observed  between  front  and  back  of  the  sample. 

5.  Model  performance  compared  to  experimental  results 

Entropic  heat  generation  terms  were  measured  as  a  function  of 
the  DOD  for  each  electrode  material  as  described  in  Section  4.  These 
results  were  used  in  the  model  described  in  Section  2.  We  found  that 
heat  loss  along  the  tab  area  has  a  significant  effect  on  the  simulation 
results  for  prismatic  batteries.  The  battery  simulation  consisted  of 
solving  for  two  primary  unknowns:  the  heat  transfer  (1)  from  the 
leads/clamps  and  (2)  from  the  convective  cooling  of  the  battery 
faces.  Two  sets  of  experiments  at  the  4C  discharge  rate  were  used 
for  calibrating  these  heat-transfer  conditions:  one  with  the  bat¬ 
tery  insulated,  and  one  exposed  to  ambient  (natural  ventilation). 
The  insulated  case  was  used  to  eliminate  convective  cooling  so  that 
heat  transfer  from  the  tabs/clamps  could  be  measured.  Once  heat 
transfer  from  the  clamps  was  quantified,  the  air-exposed  experi¬ 
mental  data  were  used  to  calibrate  the  convective  cooling  on  the 
large  faces  of  the  battery. 

The  heat  transfer  from  the  leads/clamps  was  measured  by 
attaching  three  thermocouples  to  each  tab/clamp.  Thermocouples 
were  attached  to  the  clamps  at  0.5,  1.0,  and  1.5  in.  (1.27,  2.54, 
and  3.81  cm)  distances  from  the  electrode  tab/clamp  interface. 
During  the  4C  discharge  rate  experiments,  we  found  that  the  Al 
battery  tab/clamp  configuration  exhibited  a  temperature  gradient 
that  would  conduct  heat  away  from  the  battery.  However,  the  cop¬ 
per  tab/clamp  assembly  exhibited  a  reversed  temperature  gradient 
that  would  conduct  heat  into  the  battery.  The  battery  discharge 
was  terminated  at  DOD  =  0.55  (t~450s)  because  of  the  large  IR 
drop  at  the  4C  rate.  At  DOD  =  0.55,  the  temperature  gradients  along 
the  Cu  tab  were  measured  as  -10°C/2.5cm  and  +12°C/2.5cm  for 
the  Al  tab.  The  temperature  gradients  were  reduced  to  essentially 
zero  at  DOD  =  0.0.  As  a  result,  as  discharge  progressed,  there  was 
extra  heat  flowing  through  the  top  edge  of  the  battery.  Fig.  13 
compares  measured  sample  temperatures  with  simulated  results. 
Nonlinearities  in  tab/clamp  temperature  gradients  began  at  about 
300s.  The  convective  coefficients  were  C=  0.175  Won-2  °C_1  and 
n  =  0.25  (see  Eq.  (8)).  The  battery  was  cooled  after  discharge,  and 


sample  temperatures  were  continuously  recorded.  The  lower  mea¬ 
sured  temperatures  after  300  s  can  be  attributed  to  heat  loss  in 
the  battery  lead/clamp  temperatures  as  discussed  above.  For  more 
accurate  comparisons,  heat  loss  through  both  contact  areas  and 
thermocouples  need  to  be  considered. 

6.  Conclusions 

The  entropies  in  various  cathode  and  anode  materials,  as  well 
as  complete  Li-ion  batteries,  were  investigated  using  the  ETMS.  We 
found  that  the  detail  shape  of  the  entropy  curves  strongly  depends 
on  the  manufacturer  of  the  materials  even  for  materials  with  the 
same  nominal  compositions.  LiCo02  has  a  much  larger  entropy 
change  than  LiNixCoyMnz02.  This  means  that  LiNi*CoyMnz02  is 
much  more  thermodynamically  stable  than  LiCo02.  The  tem¬ 
peratures  around  the  positive  terminal  a  prismatic  battery  are 
consistently  higher  (~5°C)  than  those  at  the  negative  terminal. 
This  is  consistent  with  those  observed  by  Pesaran  et  al.  on  cylindri¬ 
cal  cells  [16].  The  effects  of  using  battery-averaged  entropies  in  the 
simulation  tend  to  overestimate  the  predicted  temperatures.  When 
all  other  simulation  parameters  are  the  same,  the  temperature  pro¬ 
file  of  a  battery  is  about  5  °C  higher  (at  the  end  of  discharge)  when 
the  average-battery  entropy  is  used  when  compared  with  profiles 
obtained  from  simulations  in  which  individual  entropies  are  used 
for  the  anode  and  the  cathode.  We  also  found  it  is  interesting  that 
no  significant  temperature  gradient  was  observed  at  the  center  of 
the  stack  (i.e.,  inside  of  the  battery)  and  on  the  outside  surface  of 
the  battery  in  an  eight-cell  battery  discharged  at  the  1 C  rate  under 
natural  ventilation  conditions. 
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